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Abstract. We review gas-phase abundances in PNe and describe their dual utility as archives 
of original progenitor metallicity via the a elements, as well as sources of processed material 
from nucleosynthesis during the star's evolution, i.e., C, N, and s-process elements. We describe 
the analysis of PN spectra to derive abundances and discuss the discrepancies that arise from 
different choices at each step. Abundance results for the Milky Way and Magellanic Clouds 
from various groups of investigators are presented; the observational results are compared with 
theoretical predictions of AGB stellar yields. Finally, we suggest areas where more work is needed 
to improve our abilities to determine abundances in PNe. 
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i ■ 1. Introduction 

o ■ 

Chemical abundances in planetary nebulae (PNe) offer a window into the nucleosyn- 
thesis and interior processes occurring in low- to intermediate- mass stars (LIMS). Abun- 
dances in individual PNe illuminate this history star by star; combined results for a 
galaxy's ensemble of PNe provide a broader view of chemical evolution. In particular, 
t-H ■ the products of nuclear burning and convective dredge-up episodes depend on proper- 
ties of the specific PN progenitor, while the a elements, largely unaffected by LIMS 
nucleosynthesis, archive the progenitor's original abundances and can reveal galaxy-wide 
ly.^ | characteristics. In this review, we describe the observational data and calculation meth- 
£vq ■ ods used to derive PN abundances, discuss uncertainties in the process, and present recent 
abundance results, comparing them with theoretical predictions from AGB models. 

o : 

| 2. Observations 

Long-slit, moderate-dispersion optical spectroscopy has been the mainstay of PN abun- 
• »"H . dance determinations for decades. The resulting 1-D spectra - sometimes divided into 
sub-spectra along the slit to yield a modicum of spatial resolution - have provided abun- 
dance data for hundreds of PNe in the Milky Way (MW) and Magellanic Clouds (e.g., 
see references to Table 1). The last decade has seen the rise of integral field spectroscopy, 
in which each pixel yields a spectrum, rather than a single intensity. Such multiplexed 
spectroscopy has allowed 2-D mapping of nebular properties (e.g., Tsamis et al. 2008f 
and provided data that can be used to construct and test theoretical PN models. 

Beyond the optical region, Spitzer has enabled us to detect infrared lines from impor- 
tant ions (e.g., [S IV] , [Ar II], and [O IV]) that are otherwise invisible; Hubble has done 
the same in the ultraviolet (e.g., N III], [N V], and C III]). Expanding the wavelength 
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regime of observed ions reduces dependence on ionization correction factors, (ICFs), 
yielding more accurate abundances, discussed further in §3, below. 



3. Data Analysis 

In the absence of a customized photoionization model, abundances in a PN are typically 
calculated empirically, using a multi- level atom program such as nebular ([Shaw fc Dufour~ 1995h 
or ELS A (Johnson ct al . 2006p . While space does not permit a full examination of all 
relevant issues in data analysis, we briefly discuss the three main steps, each of which 
presents opportunities for different choices that will affect the results. 

3.1. The Empirical Method 

After removal of the instrumental signature from the observed fluxes (itself a possible 
source of systematic error), a wavelength-dependent correction for interstellar extinction 
is applied. There are many versions of this function; among those currently in use are 
Clayton et al. (1989)[|Fitzpatrick (1999)[ and |Savage fc Mathis (1979)| De-reddened line 



flux ratios will differ systematically depending on which function is chosen. 

Transforming de-reddened flux ratios into ionic abundances requires knowledge of 
atomic data, including recombination rates for permitted lines and A-values and col- 



lision strengths for collisionally-excited lines. Data cited by Mendoza (1983) have been 
widely used over the past several decades; Liu (this volume) describes some more recent 
calculations. Use of different atomic data can produce large discrepancies in the derived 
abundances (see §4.2). An important issue that needs to be addressed as we go forward 
is that when new atomic data are published, it is not easy for the end-user to evaluate 
their quality and validity; also, "the newest data are not always the best" (see Luridiana 
& Garcfa-Rojas, this volume). 

To convert ionic abundances into total elemental abundances, one must account for 
unseen ionization stages. Without a tailored model, one must use ionization correction 
factors - the ICF is the ratio of the total elemental abundance to the sum of the detected 
ionic abundances. Most ICFs exploit coincidences in ionization potentials with ions of 
oxygen and helium, the most thoroughly-observed elements; in some cases (e.g., S), ICFs 



come from modeling results (see |Kwitter fc Henry 200~lj ). Kingsburgh & Barlow (1994) 
published a set of ICFs that remains in wide use. Despite the ad hoc nature of ICFs, 
they work well in most circumstances. However, there are instances where they do not, 
as for Ne under low-excitation conditions, leading to abundances that are too low (see 
Peimbcr t et al. 1992( Henry et al. 2011 priv. comm.). 

The bottom line is that at each step outlined above, the particular choice made will 
affect the final abundances, and can potentially lead to significant disagreements com- 
pared with results by others who have made different choices, even if they started with 
identical measured fluxes. It is therefore desirable for those who work in PN abun- 
dances to make their choices public - and perhaps even to discuss them collaboratively 

in order to disperse the chaff unrelated to true abundance behavior in PNe. 

3.2. The Abundance Discrepancy - still unresolved 
When it became possible to measure weak recombination lines from light elements, it was 
discovered that they implied abundances in excess of those derived from collisionally- 
excited lines (e.g., [Tsamis ct al. 2003), sometimes by large factors. Proposed explana- 
tions invoke chemical inhomogeneities in the form of cold, H-deficient inclusions (e.g., 
ILiu et al. 2006[) or temperature and density variations in a chemically homogeneous neb- 
ula (e.g.. ITorre s-Pcimbcrt & Pcimbe rt 2003|) . Challenges to the resolution of this issue 
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have included the lack of atomic data relevant to the low temperatures assumed for 
these inclusions (though this is being addressed - see Liu, this volume), and the absence 
of definitive observational data. Since the proposed inclusions would radiate strongly 
in infrared fine-structure lines, the spatially-resolved SOFIA observations of NGC 7009 
described by Rubin (this volume) may prove useful. 



4. Light Element Abundances in the Milky Way, LMC, and SMC 

4.1. Milky Way Disk, Bulge and Halo 

In Table 1 we present a compilation of recent abundance results for He, C, N, O, Ne, 

5, and Ar in the MW disk, bulge and halo, and for the Magellanic Clouds. All the MW 
disk abundances agree within the uncertainties. The a element abundances are close to 
solar, while He, N and C are enhanced. S appears underabundant; this sulfur anomaly, 
has long been a problem to understand (Henry et al. 2004, Henry et al. (this volume)). 

Bulge PNe exhibit He, O, Ne, and Ar abundances similar to the disk, while C and N 
appear slightly less abundant, and the S anomaly remains. The halo appears deficient 
compared with the disk in all elements measured, except for C. It should be noted that the 
halo PN sample is small and diverse; the median (or average) abundances are therefore 
not as robust an indicator of overall behavior as they are for other samples. 



Table 1. Abundance Comparisons (median 12+log(X/H)) 



PN Sample 


He 


C 


N 


O 


Ne 


S 


At 


Ref. 


MW disk 


11.09 


8.87 


8.33 


8.61 


8.02 


6.69 


6.40 


1 




11.04 




8.13 


8.61 


8.00 






2 




11.03 


8.64 


8.18 


8.62 


8.13 


6.87 


6.43 


3 


MW Bulge 


11.05 


8.44 


8.10 


8.71 


8.03 


6.99 


6.26 


4 




11.11 




8.11 


8.57 


7.93 


6.79 


6.34 


5 




11.04 




8.01 


8.48 


7.79 


6.59 


6.37 


6 


MW Halo 


10.97 


9.00 


7.54 


7.83 


6.82 


5.39 


4.78 


7 




10.99 




7.76 


8.46 


7.81 






2 


LMC 


11.00 


8.17 


7.75 


8.24 


7.62 


7.08 


5.90 


5; 8 for C 




11.01 


8.56 


7.79 


8.31 


7.56 


7.27 


6.00 


9 












7.75 


6.47 




10 


SMC 


11.11 


8.25 


7.17 


8.04 


7.40 


6.11 


5.57 


11; 12 for C 




10.95 




7.26 


8.10 


7.35 


6.72 


5.63 


13 












7.45 


6.04 




10 


Orion 


10.99 


8.52 


7.73 


8.73 


8.05 


7.22 


6.62 


14 


Sun 


10.93 


8.43 


7.83 


8.69 


7.93 


7.12 


6.40 


15 



iHenry et al. (2010)| 2: [S tanghcllini & Haywood (2 010) | 3: jPottasch fc Bernard-galas (2010" 
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4.2. The Magellanic Clouds 

As has long been known, the Magellanic Clouds are less chemically evolved than the 
MW. Both Clouds are deficient in He compared with the MW dis k, the LMC by ~0. 05 
dex, and the SMC by ~0.1 dex. The SMC He abundance from |Shaw et al. (2010)] is 
0.15 dex higher than that of Idiart et al. (2007), arising from the use of a new He + 
recombination rate (Porter ct al. 2007). This is a perfect example of why it is imperative 
to understand how results are obtained by different groups before forming conclusions 
about any discrepancies. The median N abundances in the LMC and SMC are below 
the MW disk by -0.5 dex and -1 dex, respectively. The LMC and SMC median C 
abundances are similar to each other, presumably due to the lack of hot-bottom burning 
in the SMC PNe progenitor stars. The abundances of the a elements are lower in the 
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LMC than in the MW disk by ~0.4 dex, and in the SMC by ~0.6 dex, except for Ar, 
which is lower in the SMC by ^0.8 dex (though shortcomings in the ICF for argon may 
be a contributing factor). 



5. Abundances of s-process Elements and Iron in Milky Way PNe 

The s-process elements are produced in the interpulse phase of the TP-AGB in the 
region between the hydrogen- and helium-burning shells; third dredge-up (TDU) then 
carries them, along with 12 C, to the surface. Since the hrst report of s-process elements 

much progress has been made, especially in the 
detected lines of Br, Kr, Xe, Rb, Ba, and Pb 



in PNe by Pequignot fc Baluteau (1994) 



last several years. Sharpee et al. (2007) 



and possible lines of Te and I, and found that abundances of Kr and Xe were enhanced. 
Sterling fc Dinerstein (2008)| ;SD08) detected Kr and/or Se in 81 of 120 PNe; in non- 



Type I PNe they found average enhancements of ~0.3 dex for Se, and ~1 dex for Kr. Type 
I PNe exhibited little s-process enrichment, which SD08 attribute to smaller intershell 
masses and efficient dilution of s-process material in the massive AGB envelopes. They 
found a positive correlation between s-process enrichments and C/O, as expected in 
TDU. Future progress in learning more from s-process abundances in PNe relies heavily 
on knowing transition rates and collision strengths (See Sterling, this volume). 



Delgado Inglada et al. (2009) observed [Fe III] lines in 33 low-ionization PNe. They 



found a median 12+log(Fe/H) = 5.85, and a range in depletion, [Fe/H], from -1.01 to 
-3.2, implying depletion > 90%, and Md ust /M gas ^ 1.3 x 10~ 3 . Further, they found that 
the depletion correlates somewhat with C/O ratio (see Delgado Inglada, this volume). 



6. Abundance Correlations and Model Tracks 

6.1. Helium, Carbon, Nitrogen & Oxygen 

Fig. 1 shows N vs. He for a variety of PN samples. The MW disk PNe are taken from 
Henry et al. (200*4)] reanalyzed with consistent atomic data; these PNe all have good 



measurements of T\pni)i Trjy //], and N e rgjn, so the abundances are well determined. 
Type I PNe in the MW and LMC occupy the upper right of the diagram. Note that the 
plotted sample of SMC PNe lacks any Type I PNe. In this and following plots we have 



overlaid predicted PN abundances taken from Karakas (2010) and Marigo (2001) for the 
metallicities of the MW, LMC, and SMC, for initial masses from 1-6M and 0.8-5M Q , 
respectively; solar values are shown by vertical and horizontal lines. The important point 
to note is that the models span the actual observations. 

Fig. 2 plots C/O vs. O for various objects. The fairly narrow trend outlined by H II 
regions, damped Lyman-a systems, halo stars, and disk stars is completely ignored by 
PNe, their wide swath testifying to their non-uniform C production. 

In Fig. 3 we see N/O vs. O/H for the same samples as in Fig. 1. Note the trend of 
decreasing N/O with increasing O/H exhibited by the LMC PNe, and also by the MW 
Type I PNe; this is a sign of hot-bottom burning at the base of the convective envelope 
in massive TP-AGB stars, converting C and O into N. Note also that neither the SMC 
nor the MW non-Type I PNe exhibit this tendency. 

6.2. Neon 

In Fig. 4 the Ne and O abundances are plotted for H II regions and blue compact 
galaxies (H2BCG; Milingo et al. 2010) in addition to PNe in the MW disk and in M31's 
disk (jKwitter et al. 2011)) . The tight relation evident for the first two is followed, by PNe, 



but with more scatter. Model predictions from Karakas (2010) are shown as well. The 
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Figure 1. N vs. He with model predictions. Type I PNe occupy the top right of the diagram. 







i 


o 




Halo Stars (Akarman etal. 2004) 
MWG Disk PNa IHBnry etal. 2004) 
A DLfe(Cookeet a l.201l) 
□ LMC PNe (Slannhellini el al. 2005) 


o o 

° o 
° o °°° 


o 




- 


k 


- a ^ - ,411*- 

A A^A* / .» ♦^•••» 

♦ v.»vs ♦ ♦ L 

□ c 

1,1 



Figure 2. C vs. O. Solar values are indicated by dashed lines. 
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Figure 3. N/O vs. O/H with solar values indicated by solid lines. Note the tendency in MW 
Type I and LMC PNe of decreasing O as N/O rises, a result of hot-bottom burning. 

tendency for MW PNe to lie slightly, but systematically, above the H2BCG line could be 
interpreted as an indication that PNe are more enriched in Ne than H II regions of the 
same metallicity (but see §3.1). Also, Fig. 5 shows Ne/O vs. O for the same PNe samples 
as in Fig. 1 and there is no compelling evidence for a trend, which indicates either that O 
is not altered in LIMS nucleosynthesis at MW or MC metallicities, or that both Ne and 
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O are changed similarly; Maciel et al. (2010) find the same result. Wang & Liu (2007) 
concluded that Ne and O production only occurs for 12+log(0/H)<8, below even the 
SMC, indicating that the first possibility is likely the correct one. 



H2BCQ 

LSFto H2BCG 

AQB Models (Karakas 2010) 

MW Disk PHe (Henry et al. 2004) 

SMC (Shaw el al. 2010) 

LMC (Leisy & Dennifeld 2006) 

M31 Disk (Kwitter el al. 2011) 




Figure 4. Ne vs. O, with models. Solar values are indicated by dashed lines. For H2BCG 

references, see Milingo et al. (2010) 
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A LMC (Chiappini et al. {2009) 
■ SMC (Shaw etal. 2010) 



-2.000 

7.000 7.500 



12+log(0/H) 

Figure 5. Ne/O vs. O/H. The lack of a discernible trend implies that either O is not affected 
by LIMS evolution, or that both Ne and O are similarly altered. 



6.3. Sulfur and Argon 

In Fig. 6 we plot S vs. O, with other objects and model predictions as in Fig. 4. Also 
shown are the predictions (marked by X) of photoionization models for five MW PNe, 
confirming that they exhibit the sulfur anomaly. There is a new possibility as to how the 
sulfur anomaly might be resolved, having to do with the existence of a S +4 zone predicted 
by RHD models (see Jacob, this volume); we will be eager to see how it succeeds. 

Fig. 7 shows Ar vs. O plotted for the same objects as in Fig. 4 for Ne. The trend set 
by the H2BCG sample is well followed by PNe, though the scatter is larger. The general 
concordance of Ne and Ar with O bolsters the view that these a elements do grow in 
lockstep, as predicted by nucleosynthesis arguments. 



6.4. Morphological Correlations 

Striking correlations emerge when comparing different populations of PNe. Maciel & Costa (2010) 
show in their Fig. 1 that N and He abundances are significantly higher in asymmetric 
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Figure 6. S vs. O, with models. Solar values are indicated by dashed lines. 
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Figure 7. Ar vs. O. Solar values are indicated by dashed lines. 



PNe compared with symmetric PNe; the opposite is true for C, implying that asymmetric 
PNe originate from systematically more massive stars. In addition, O, Ne, Ar, and S are 
all slightly higher in asymmetric PNe. This behavior is to be expected if the higher-mass 
progenitors of asymmetric PNe are younger and were born from a more-enriched ISM. 



For the MCs Stanghellini et al. (2009) found that in their samples, the C abundances in 



SMC PNe are significantly higher than in LMC PNe (their Fig. 13). 



7. Conclusions 

We end with a look ahead. Atomic data calculations for low temperatures and new 
elements and ions are now being carried out. A small group of investigators has begun 
a systematic look at the empirical abundance process as described in §3.1, with an eye 
toward evaluating and possibly even standardizing some of the procedure. Studies on the 
dependence of ICFs on metallicity, geometry and nebular ionization are also proceeding. 
Ongoing 2-D nebular abundance mapping and kinematic studies will enable 3-D nebular 
modeling, with the potential to solve, or at least to shed light on the abundance discrep- 
ancy. We anticipate exciting and illuminating new results at the next symposium. 
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